TOWARD THE FOURIER LAW FOR A WEAKLY 
INTERACTING ANHARMONIC CRYSTAL 



CARLANGELO LIVERANI AND STEFANO OLLA 

Abstract. For a system of weakly interacting anharmonic oscillators, 
perturbed by an energy preserving stochastic dynamics, we prove an 
autonomous (stochastic) evolution for the energies at large time scale 
(with respect to the coupling parameter). It turn out that this macro- 
scopic evolution is given by the so called conservative (non-gradient) 
Ginzburg-Landau system of stochastic differential equations. The proof 
exploits hypocoercivity and hypoellipticity properties of the uncoupled 
dynamics. 



1. Introduction 

The problem of deriving heat equation and Fourier's law for the macro- 
scopic evolution of the energy from a microscopic dynamics of interact- 
ing atoms (hamiltonian or quantum), is one of the major goals of non- 
equilibrium statistical mechanics [2]. 

Although we are still very far from a rigorous mathematical derivation, 
we have now some understanding of the needed ingredients. 

It is clear that heat equation is a macroscopic phenomenon, emerging af- 
ter a diffusive rescaling of space and time. It is also clear that non-linearities 
of the microscopic dynamics are necessary, since in a linear system of inter- 
acting oscillators energy may disperse ballistically and thermal conductivity 
results infinite |10| . Non linearities of the interaction should give enough 
chaoticity and time mixing such that locally the system, in the macroscopic 
time scale, is in a state of local equilibrium. This should be intended in 
terms of a scale parameter e: in a region of linear size e, at a large time 
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scale e~ b t, the system should be close to equilibrium with temperature given 
by the local average of kinetic energy. This statement of local equilibrium 
should be intended in the stronger sense that locally the dynamics is close to 
an equilibrium dynamics. Since energy is a conserved quantity, it can only 
evolve by moving between different regions of linear size e~ l through energy 
currents. Because of the size of the regions and the fact that in equilibrium 
energy currents have null average, one should look at time of the order s~H 
in order to see some exchange of energy between boxes at different tem- 
perature. In other words a central limit theorem for the energy currents 
is involved, and the thermal conductivity is then given by the space-time 
integral of the current-current correlation (Green-Kubo formula). This con- 
ductivity will be convergent if the system in equilibrium has enough mixing 
properties. 

To perform the above program, in a mathematical rigorous way, from 
a purely deterministic Hamiltonian dynamics, it is at the moment a too 
difficult challenge. 

In the last years some mathematical results have been obtained by per- 
turbing the dynamics with energy conserving stochastic forces. The purpose 
of this stochastic perturbations is to give the ergodic and chaotic properties 
to the system without modifying the macroscopic behavior of the evolution 
of the energy. 

This strategy has proven successful for systems in the hyperbolic scaling 
(b = 1), when momentum conservation is also preserved by the stochastic 
perturbations [9J, obtaining Euler system of equations for compressible gas 
as macroscopic equation, at least in the smooth regime. 

In the diffusive scaling (b = 2) this problem is still very challenging even 
in presence of the stochastic perturbations. The difficulty is essentially 
involved in the space-time rescaling and the corresponding central limit 
theorem. 

In this paper we develop a weak-coupling approach to the problem of 
energy diffusion that permits to separate the time limit from the space one. 
We consider a finite system of anharmonic oscillators, whose hamiltonian 
dynamics is perturbed by a noise that conserves the kinetic energy of each 
oscillators (we consider oscillators that have at least 2 degree of freedom). 
The noise could be though as modeling some chaotic internal degree of 
freedom of each atom. 

The oscillators are weakly coupled with a small parameter e. Conse- 
quently the exchange of energy between oscillators is given by the currents 
associated to the hamiltonian mechanism, but multiplied by e. The noise 
drives each atom towards (microcanonical) equilibrium where currents have 
null average. In a time scale of order e~ 2 current fluctuations are able to 
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move energy around the lattice if gradients of energy are present between 
atoms. 

We prove in fact that in the limit as e -*■ 0, in the time scale e~ 2 t, the 
energies of the atoms evolve autonomously following the solution of a sys- 
tem of stochastic differential equations, conservative of the total energy (cf. 
( 13. 8p ). It turns out that this macroscopic stochastic evolution has already 
been considered in the hydrodynamic limit literature and it is called non- 
gradient Ginzburg-Landau model [12]. Consequently, using the techniques 
developed in [12], one could try to prove that under space-time diffusive 
rescaling, the energy evolves following a non-linear heat equation. The 
results in [T5] do not apply directly to ( 13. 8 p due to degeneracy of the coef- 
ficients, but hopefully they can be adapted to the present situation, for the 
moment we postpone this problem. We have thus reduced the derivation 
of the heat equation to a two step procedure of which this paper rigorously 
accomplish the first step. 

We should remark here that an extension of the non-gradient technique 
of [12] directly to our original microscopic stochastic dynamics would be 
much more challenging, as this dynamics is very degenerate. 

The reason of the name non-gradient comes from the fact that the cur- 
rents of the macroscopic dynamics are not gradient, i.e. are not a given by 
the spacial gradient of a local function of the configurations of energies. It 
is interesting to note that in the purely harmonic case, the macroscopic dy- 
namics (13. 8p become gradient (see appendix |E|) . This implies a connection 
between the non-gradient property and non-linearity of the microscopic dy- 
namics. We also notice that in the purely harmonic case, because the pres- 
ence of the energy conserving noise, the microscopic energy current have an 
exact fluctuation-dissipation decomposition in a gradient plus a fast fluctu- 
ating term (see formula (lE.ll in appendix[E|), that has been already exploited 
in [T] to obtain Fourier's law. 

The main strategy of the proof is similar to other averaging principles 
([HOI?]): at a large time scale the dynamics of each atom is close to 
the equilibrium dynamics parametrized by its own energy. Energies of the 
atoms are our slow variables and evolves through their currents. But a 
simple averaging of these currents (that would occur in a time scale of order 
s^ 1 ) would not move any energy, since currents have null averages respect 
to all equilibrium measures (microcanonical). This forces us to look at 
the time scale e~ 2 , when the energy evolves due to the fluctuations of the 
currents in equilibrium. Thus we must establish a central limit theorem for 
the energy currents in the uncoupled dynamics, i.e. we have to study the 
Poisson equation 
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where Lq is the generator of the uncoupled dynamics, and j is the energy 
current between two particles. In order to prove our theorem, we need 
existence and regularity of the solution u of this equation. The generator 
Lq turns out to be hypoelliptic on each microcanonical energy surface, that 
provides regularity on the tangent direction of this surface. Yet, as the 
energy is exchanged from one particle to the other, we also needed to prove 
regularity in non tangential directions. 

For the existence of u, we prove a spectral gap in a proper Sobolev space, 
with an adaptation of hypocoercivity techniques [IB] . These techniques 
provide a precise control of this spectral gap with respect to the energy, a 
control especially needed at low energies. Such detailed informations are 
necessary in order to perform the closure of the macroscopic equations. 

The content of the paper is as follows. 



Let us consider a region A c Z d , set N = |A|, the number of sites in 
A. At each site we have a v- dimensional, v > 2, nonlinear oscillator and 
we assume that such oscillators interact weakly via a non-linear potential. 
Such a situation is described by the following Hamiltonian in the variables 
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2. The model 
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ieA Z ieA z |i-j[=l 

where € C°°(R. V , E). We use the convention £|i-j|=i = EieA E{jeA: [i-j|=i}- 
For simplicity we assume that U(q) = U(\q\ 2 ), U e C°°(R,E), U(0) = and 

cr 1 < J7 < c for some finite positive constant c, in particular this implies that 

U is radially symmetric and strictly convexo Also, we assume ||W(g)|| 2 < 

cU(q) and V(-q) = V{q). 

For simplicity of notations, we choose v = 2. All result stated in this 

paper are valid for general v > 2, with slight modifications of notations. 
In addition to the Hamiltonian dynamics, we consider random forces that 

conserve the single sites kinetic energies, given by independent diffusions on 

the spheres ||p;|| 2 = cost. In order to define such diffusions, consider the 

vector fields 

Xi ■= p}d p 2 - pid p i =: Jpi -d Pi , J = 
and the second order operator 

ieA 

The generator of the process we are interested in is then given by 

L eA =:A £ + ( r 2 S (2.1) 

where A e = {H^,-}, is the usual Hamiltonian operator and a > measures 
the strength of the noise. Clearly, L £j a is the generator of a contraction 
semigroup P\ in L 2 (M. 2N ", m e ) with stationary measure m e and P^H^ = H^, 
for all t e M + . Next, we must specify the initial conditions. 
The Gibbs measures at temperature (5^ are defined by 

mP(dq,dq) = Z^)e^ H ^dpdq 

and are the stationary (equilibrium) probability measures for the dynamics 
(the canonical ensemble) . As reference measure we pick the one correspond- 
ing to ft = 1 and we denote it by m £ . Notice that for e small enough, m £ 
and the product measure m are equivalent. To simplify notations, we also 
assume that U is such that 

ZJl) = 1 = f e-^ 2 / 2+U ^dpdq. 
JR 4 

1 The general, non radial, case can be treated exactly in the same way at the only price 
of a much messier algebra. On the contrary, the non convex case could hold interesting 
surprises and hopefully will be investigated in the future. 
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We assume that the system is started in an initial distribution du := 
F £ dm £ = F dm . 

Assumption 1. We assume that Fq e L 2 (M 4iV , mo). 

For each T > 0, the Markov process just described defines a probability on 
Q = C°([0,T],M 2 H A l). We will use co t = (q(t),p(t)) to designate the elements 
of f2 at time t. 

Remark 2.1. In the following we will suppress the subscripts and super- 
scripts A, when this does not create confusion. 

Remark 2.2. Here we have free boundary conditions. It should be possi- 
ble to treat more general stochastic boundary conditions (e.g. having the 
particle at the boundary perform an Ornstein- Uhlenbeck process at a given 
temperature) by the same method, we avoid such a generalization to simplify 
the presentation. 



3. The results 
The single particles energies are 

£f(q,P) = T:\Pif + U{q,i) + -e Y, V(Qi-Qj)- 
A 1 |i-j|=l 

The time evolution of these energies is given by: 

d£f 



e 



\i-k\=i 

where the energy currents are defined by 

ji,k = ^V(qi-q k )-(pi+p k ). (3.2) 

Note that ji,k = ~jk,i an d that they are functions of the qi,Pi,q k ,p k only. 

If e = the dynamics is given by non-interacting oscillators, and con- 
sequently the energy of each oscillator is a conserved quantity. So for 
e = there is a family of equilibrium measure parametrized by the vec- 
tor a = (aj)j e A of the energy of each oscillator. This is given by the 
microcanonical measure associated to the Hamiltonian flow Hq on the sur- 
face 

£„:={g,p = a l =S?(q,p) = ^\\p t \\ 2 + U(q t )}=X teA Z ai . (3.3) 

Clearly, letting /i a be the microcanonical measure on the 3 dimensional 
surface E , we have /i£ = ®ieA/J- ai - By the symmetry between p and -p it 
follows that fig,(ji,k) = for each a. 
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We are interested in the random variables determined the time rescaled 
energies 

£f(t)=Ef(q{e-H),p{e-H)). (3.4) 

In order to define the parameters of the mesoscopic evolution, consider 
the dynamics of 2 non-interacting oscillators (e = 0), each starting with 
the microcanonical distribution with corresponding energy a% and a 2 . Let 
us denote by E ai)02 (-) the corresponding expectation in this equilibrium 
measure. We will show that the following function on M.+ 

E ai , a2 0i, 2 (t)ji, 2 (0)) dt, (3.5) 

o 

is well defined. More, in Lemma 18.21 we prove that 

7 2 (ai, a 2 ) = a^G^di, a 2 ) (3-6) 

where G is a positive symmetric smooth function. Correspondingly we 
define the mesoscopic current by the antisymmetric function 

a(a 1 ,a 2 ) = e u ^(d ai -d a2 )(e- u ^ 1 2 (a 1 ,a 2 )). (3.7) 

where U(a) = - £ ■ log£(aj), and Z(a) is the energy density distribution 
under mo, that behaves like a for small a. 
Here is our main result: 

Theorem 1. In the limit e -> ; the law of {£ fjjeA converges to the weak 
solution of the stochastic differential equations 

d£i= Y, a(£i,£ k )dt+ £ j(£i,£ k )dB {i;k} (3.8) 

k'\i— k\=l fc:|i— fc|=l 

with = _ -B{fc,i} independent standard Brownian motions. Where the 

law of £i(0) is given by the marginal of F e dm £ on the £ £ . 

Notice that the generator of the diffusion (I3.8P on is given by 

£= Y, {l(£ l ,£k) 2 (d £i -d £k r + a(£ l ,£ k )(d £i -d £k )). (3.9) 

|fc-i|=l 

Since 7 2 and a are smooth function on , the uniqueness of the weak 
solution of (I3.8p follows by applying the results in [3]. 

Observe that there is a family of product probability measures 

Y[Z(a i )e-^N(/3)- 1 da i = Y\e- {l3a ' +uia ' )) N(/3)- 1 da u /3 > (3.10) 

ieA ieA 

that are stationary and reversible for the diffusion generated by (13.91) . 

As we will see shortly the proof of Theorem [1] relays heavily on the fact 
that the unperturbed microscopic dynamics P l K generated by L ,a has strong 
mixing properties (hypocoercivity) . This is itself a non trivial result which 
we believe worth stating separately. Let %\ be the Sobolev space of order 
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one on Eq with respect to the microcanonical measure /x„ and a properly 
rescaled Riemannian structure (see Section [6] for more details). 

Theorem 2. For each set A c Z d ; there exists C,r > such that, for each 
energies a e (0, oo) A and o e (0, 1), i/ie following holds true 

• The semigroup Pq ^ contractive in L 2 (Xq,, /%)■ 

• For eac/i smooth function f e H„ suc/i t/iat fJ>a(f) = 0, Zioids 

linVII^^^II/ll^- 

The proof is given in section [61 

Before discussing the proof of the above results let us indulge in several 
remarks. 

Remark 3.1. By (13.71) we can rewrite the generator as 

C= £ deJeVM-fieMfa-deJ (3.11) 

|fc-i|=i 

Remark 3.2. The process (13. 8p is close the the one studied by Varadhan in 
[12\, yet it is not covered by such results (due to the degeneracy at zero of 
the diffusion coefficients and the non strict convexity of the potential of the 
invariant measure). In any case, the extension of Varadhan 's work to the 
present case would allow to obtain the heat equation in the present setting 
via a diffusive scaling limit of space and time fhydrodynamic limitj. 

Remark 3.3. Note that both 7 and a depend on a. One can wonder if 
equation (13.81) does admit a limit for small noise. Indeed, if U and V are 
quadratic then both 7 2 and a are proportional to a~ 2 , see (1E.4j) . Hence 
the energy exchange for small noise is faster than the time scale we are 
exploring. This is due to the fact that in the quadratic case the solutions 
are quasiperiodic. On the other hand, for each U,V a positive measure of 
such quasiperiodic solutions will survive by KAM theorem at least for small 
energies so it may be possible that a small noise limit exists, upon rescaling 
time, as in the quadratic case. The present results allow only upper bounds 
(which agree with the quadratic case), but it is unclear if a sufficiently exact 
scaling still exists. 

4. Proof of Main theorem 

This section is devoted to proving Theorem [T] by using several results 
detailed in the later sections (more precisely we assume Proposition 16.11 and 
Lemmata 17. 3[ 17.41 18.11 and 18.21) . Our strategy is the first establish tightness 
and then to show that the accumulation points satisfy (13.81) . Since (13.81) has 
a unique solution, the process have a unique accumulation point whereby 
proving the existence of the limit. 
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4.1. Tightness. Here we start studying the processes {Sf(t)} defined in 

(EH). 

Lemma 4.1. There exists £q > such that for each T > 0, the processes 
{£f(t)} t <T, < e < So, are tight. 

Proof. The proof use a standard backward/forward martingale decomposi- 
tion argument (cf. [TTj IS]). We recall it here. 

Let us start the process with the equilibrium distribution m e . Then 
the time reversed process, in a given time interval, is a Markov process 
with generator given by the adjoint L* = -A e + a 2 S. Remark that, since 
X^p k = -p k , we have Sji )k = -ji,k- So we can decompose 

Jets' 2 £ rts~ 2 £ rte~ 2 

ji,k( S ) ds = / L eji,k( S ) ds ~7r^ L *eJi,k( S ) ds 

o 2a z Jo la* Jo 



2a 2 t£ 2a 2 t£ ~ 
where are continuous martingales, adapted respectively to the forward 
and backward filtration, that can be represented by the stochastic integrals 

Mt = f\x i3j , k )(s)dwf(s) + J\x k j hk ){s)dwt{s) 

where w^(t) and w^{t) are standard Wiener processes adapted respectively 
to the forward and the backward filtration. Consequently the tightness 
follows from the tightness of each of these stochastic integrals. Noticing 
that X k jj^ k = Jp k ■ VV(qj - q k )/2 is in L p (m e ) for any p < oo, by Doob's 
inequality: 

E me (sup (Mt ) 4 ) < (^)V E ((M^) 2 ) 2 < CT 2 
This imply that the Kolmogorov criterion for tightness is satisfied: 

supE m£ (|^( S )-^(t)| 4 )<^(t-,) 2 

e a 

In non-equilibrium, because of the assumption that the initial distribution 
F is in L 2 (m £ ), the above argument extends immediately by a simple use 
of Schwarz inequality. □ 

Observe that with the same argument we can also establish the following 
bound for any p > 1 : 

- rt ~\P 

/ ji,k(s) ds 
.Jo 



E me sup 

\0<t<T 



CTp/ 2 



where C is a constant independent of £,T,o~. 

Once we have the tightness all is left to prove is that the limit is unique, 
this is the content of the next section. 
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4.2. Identification of the limit. The goal of this section is to prove that 
any accumulation point of the laws of {£f (t)}^\ must satisfy equation (13.81) . 
More precisely, since 



Sf{t)-Sm= £ e f * 

fc:|fc-i|=l J0 



(4.2) 



we want to show that, for every couple i, k e A such that \i - k\ = 1 there 
exist orthogonal martingales -M £ ik = ~<Mj }i with zero average and quadratic 
variance given by 

2 [\ 2 (£?(s),£ £ k (s))ds (4.3) 
Jo 

and such that, for each t > 0, 



limEl 

£^0 1 



re~ A t r~t 

£J Q ji,k( s ) ds -J Q a{£!{s),£i{s))ds + Ml k (t) 



0. (4.4) 



To prove (14.41) it turns out to be useful to introduce a cutoff for high 
energy = Xk(£°), where \k '■ -*■ [0, 1] are smooth positive functions 
with support in [0, K + 1] A , such that Xk(o) = 1 if o e [0, i^] A . Then define 

By the symmetry between p and -p it follows that = for each 

a e M A and any K > 0. Arguing as in estimate (14. li . we have that 

n 2\ 



lim supE m£ I sup 



K^oo £>0 



, 0<t<T 



\ JiA 1 ~€k)(s) ds 
Jo 







By Proposition I6.1[ it follows that the equation 

L^Ui^k = ji,k 



(4.5) 



(4.6) 



has a unique solution with zero average with respect to all measures /v 
Note that - -jk,i, thus also u^k = -Hfc.i- 

In addition, Lemma O implies u hk e C°°([M 4 \ {0}] 2 ,M) n C°(M 8 ,M). 
Observe that, since L conserves the energies of all particles, denoting uf- = 
Uij^K we have 



T i, K - n K 



(4.7) 



We can thus write L £ = L + eL* with 



L *f=o E W( qi - qj )-(d Pi -d Pj )f= £ VV( qi - qj )-d Pi f, (U 
/ |i-j|=i \i-j\=i 



So we have 



Leufk - + eL*uf k , 



K 



(4.9) 
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Hence, denoting a path by u s = (q(s),p(s)), we have 
E Jo &( s ) ds = e J o [ L eU^ k -eL^ k ](u s )ds 
= (eu« k {uj e - H ) -eu« k (uJo) -M*£{t) - e 2 £ * L*u« k (u s ) ds^ 
where we have introduced the centered martingale 

M tf(t) ■■= ^U^e-H) - eu« k (u ) - e £ % L e u« k (u> s )ds. 
Due to the property Ui k - -u k i we have Mf' k K = -Mf', and the same for all 



the derived martingales. 



The quadratic variations of M^, K is given by 



jeA J 



(4.11) 



Xe t 
i 2 K {X j u^ k ){XjU i , )k ,){u s )ds. 

In order to close the evolution equations, as e -> 0, we need to prove that 



e^O 

and that 







(4.12) 



lirnK 

£^0 



Xe t 
He [(XjUi,k)(XjUi> >k >)M - ^^((XjUi^iXjUi,^))] ds 

(4.13) 

These are consequence of Lemma I7.4[ 17.31 and the following lemma: 

Lemma 4.2. Let f(cu) a function in L 2 (m ) such that Ha(f) = for all 
a e \ then 

1 r 1 : fMd £ 

Jo 



limE* 

s->0 r ' 



0. 



Proof. By using Jensen's inequality and stationarity, it is enough to prove 
that 

•T ri ^ 2\ 

x ts 



lim limEi, 

T^oo e^O mE 



1 r T 

T Jo f( Us ' di 







(4.14) 



Since lim e ^ E^ £ = E^ o , i.e. the expectation with respect to the dynamics of 
the non-interacting oscillators starting with the product of canonical mea- 
sures mo, which are convex combination of the microcanonical ones. Then 
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the result follows by the ergodicity of the dynamics of the single oscilla- 
tors. □ 



So far we have obtained that 



e^O 



; £ H (e-*j$,(8) - /^(L.u*)) ds - Mlf(t) 



0, (4.15) 



where the martingales M^' k K have quadratic variation given by 



{{Mlf ,M*;J»(0 = 2a 2 £ Zk(£ (s)) 2 ^(s) aX jUi , k )(X jUil , kl )) . 

je{i,k}n{i' ,k'} 

(4.16) 



Next, we remove the cutoff on the energies. Observe that 

L* u f,k = ^KL*u ijk + u ik L^ K , 
which, writing £^(a) for 9 aj ^(a), implies 

/j,a(L*uf; k ) = £K(a)Vg,( L * u i,k) + 2 £ fx (a)/^a(«i,fe VV(gj - qj) -pj) 

\i-JH 

= €K(a)l*a(L*u iik ) + 2^P(a)fia(ui ik VV(qi-q k ) -p k ) 

= ^ K {a)^a{L^ k ) +^ ) (a)7 2 ( a i) a fc) 

Since m £ (£j > K) is exponentially small in K, and 7 2 does not grow faster 
than poly normally, we have 



sup 



Xe t 
(v£°(u s )( L * u *k) -iK(£?(u s ))^ {Wa) (L*u ijk ))ds 

< Csup £m e ( 7 2 (£ K 



2 — o. 



Since the above limit and (14.51) are uniform in e, and since XjUi )k is in 
L 2 (m £ ), we can take the limit K -*■ oo in (14.151) and obtain 



limE^ 



J o 3i,k{s)ds-e 2 ^ {u]a) {L*u i:k )ds + M £ i k {t) 



0, 



(4.17) 



where M? k (t) is a martingale whose quadratic variation has the expression 
given by ( I4.16|) with £ K substituted by 1. 

Lemma 4.3. For each a k ,qi,pi setting 

uMuPi) ■= J u kii (q k ,p k ,q h pi)iJ, ak (dq k ,dp k ) 
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holds true 

ik(qi,pi) = J V(qi-q k )fi ak (dq k ,dp k )- J V{(^-q k )^a k {dq k ,dp k )^ S Q{d^dp'^. 

In particular, U{(qi, -pi) = and XiUi = 0. 

Proof. First of all note that, due to the symmetry between p and -p, 
J jk,i(Qk,Pk,qi,Pi)^a k ( d( lk,dpk) = J d qi V(qi- q k )pi/j, ak (dq k ,dp k ). 
Thus, setting V{q j ) = J V{qi - q k )fi ak (dq k ,dp k ), 

J jk,i(lk,Pk,qi,Pi)l^a k (dq k ,dp k ) = L 0j{ i } V(qi,Pi). 
In addition, by the product structure of the generator, 

jk,i = io,A«it,! = Lo,{i}U k) i + L ,{k}U k j. 

Integrating the above we have 

L ,{i}V(qi,Pi) = J [L Q){i} u kti + L 0:{k} u k:i ]fi ak (dq k ,dp k ). 

Since fi ak is the invariant measure of L ^ k j, we have 

Lo,{i}V(qi,Pi) = L 0i{i} Ui. 

By Proposition 16.11 applied with A = {i} it follows that the only solutions 
of the above equation are of the form ui = V + for some function /. 

Next, since is of zero average by construction, f(£f) = -p,go(V). □ 

Thus, if / ^ {i, k} we have 

dp i Ui ! k-W{q i - qi) Ha k (dq k ,dp k ) ^{dqi^dpi) fi ai (dqi,dpi) 

= JJ dpfii ■ VV(qi-qi) /J, ai (dqi,dpi) /2 ai (dq h dpi) = 

due to the antisymmetry of d Pi Ui with respect to Pi established in Lemma 
14.31 From this follows 

lia{L*u i)k ) = Ha\yV{qi- Qk)(d Pi -d Ph )ui tk ] =■ a(a h a k ). (4.18) 

In fact, again by the product structure, the above is a function of ai,a k only. 
It is also convenient to define 

7 2 (ai,a fc ) = cr 2 ^a({XiU i)k ) 2 + {X k u iyk ) 2 ) = -fia(u i)k ji )k ) ■ (4.19) 

Accordingly, by ( I4.16P and Lemma I4.3[ 

{(Ml k ,Ml k ,))(t)=2a 2 £ ^( s )((X jUi:k ) 2 )(6 i=ilik=k/ - <W,fc=i') (4.20) 

je{i,k} 
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Lemma 4.4. 

limE m£ p (a(£?(u E -. s ),£° k (u E - 2s )) - a(£^(s),£ e k (s))) ds 

e^O Jo 

and similarly for j 2 . 

Proof. By stationarity and Schwarz inequality we have 



E, 



<tm £ (\a(£?,£° k )-a(£f,q)\ 2 ) 



Since we prove in Lemmata 18. 11 that 01(01,02) is uniformly Lipschitz in 
K 2 the result follows by elementary arguments. □ 



Applying Lemma S3 to (14TT71) and (gjZOD we obtain that 



liniZ, 



f ji, k (s)ds- ra(£t(s),£i(s))ds + Ml k (t) 
Jo Jo 



0, (4.21) 



which, remembering (I4.19p . yields the wanted result in equilibrium. Our 
general claim (14.41) follows by Schwarz inequality. 

5. Hypoellipticity and regularity on energy shells 

We start here by studying the single site dynamics. Here p = (p 1 ,p 2 ), q = 
(o 1 ,^ 2 ) will be the coordinates. We will use L to designate the generator 
(12.11) for e = and A consisting of only one site. Since L preserves the 
energy, we can look at its action on each energy surface 

IpI 2 



S a = {(q,p) € 



+ U{q) = a}, a > 0. 



(5.1) 



Lemma 5.1. For each a > the operator L is hypoelliptic on E a . 

Proof. We must study of the Lie algebra generated by the vector fields 

Co = X = Jp-d p , 

B = A Q =p-d g -VU(q)-d p . 
We obtain 

Ci =: [C , B] = Jp ■ d q + JVU(q) ■ d p . (5.2) 

[C U B] = 2 JVE%) • d q - {D 2 U{q)J + JD 2 U{q)}p ■ d p . (5.3) 
By our assumption on U 

D 2 U{q)=AU'\q)q®q + 2U\q)t. (5.4) 

Since 

J{q®q) + {q®q)J =\q\ 2 J, 
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we have 

D 2 U(q)J + JD 2 U(q) = 4{U" (q)\q\ 2 +Tl'}J =: ((q)J. 
Finally we define C2 by the relation 

[d,B] = 4U\q)Jq ■ d q - 4{t7"(g)|g| 2 + Tl'}Jp ■ d p =: 2C 2 - ({q)C . 
Observe that, setting M = (\p\ 2 + \ VU(q)\ 2 ) 1 / 2 , the vectors 

z! = j±, z = ^, Z * = c ^jr (5 - 5) 

form an orthonormal base of the tangent space of E a , hence the Lie Algera 
generated by {Co, B} spans the space of any energy shell S a . This concludes 
the proof of the hypoellipticity of L. □ 

By the above results trivially follows the claimed hypoellipticity. 

Lemma 5.2. For each A c Z d , > ; i e A the operator L 0; a is hypoelliptic 
on S„. 

By Hormander theorem [B] hypoellipticity implies that if there exists a 
solution u for the equation 

L 0A u = g (5.6) 

where g is a C°° function when restricted to any energy shell, then also u is 
C°° when restricted to any energy shell. 

6. Hypocoercivity 

We will prove the existence of the solution of equation (14. 6 j) by proving 
a spectral gap for the generator Lo on each energy shell £„ in a proper 
Hilbert space. More precisely we consider the Hilbert spaces determined by 
the scalar products 

(h,g)jq ■= \A\{h,g)a + Y J T> aA h ^9) 

(6.1) 

T>a,i{h,g) := Y,( C k,l h > C k,l9)a + {Bih,Big)a. 

where (h,g)a = / s hgdfia and the C^i,Bi are the vector fields defined in 
section [5] relative to the particle (qi,pi). 

By a slight abuse of notations we will use Pq to denote the strongly 
continuous semigroup generated by Lq both in L 2 and for each a. 

Note that the above norm is equivalent to the standard Sobolev space Ha 
on the Riemannian surface Sq where the Riemannian structure has been 
rescaled to have the diameter of each S a . equal one independently of the Oj. 
More precisely, there exists c > such that 

cl/lk^l/U^c-^AII/lw,, (6.2) 
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Accordingly, Theorem [2] is a direct consequence of the following. 

Proposition 6.1. The semigroup Pq is contractive in L 2 . In addition, for 
each A c M. d , |A| < oo ; there exists C, r > such that for all a € (0, oo) A and 
smooth function f , such that /%(/) = ; holds 

\m\Ui<Ce- Tt \\f\Ui- 
Proof. The contractivity in L 2 follows by 

J t W P of\\l * "2 Z( C o,iHL C ,iPtf)a < 0. 

L 

To prove the second part of the theorem we use the results of Appendix 
IB] Note that for all smooth h,g,ip we can write 

m (ip(S)fj,£_(hC g)) = m (ip(S)hC g) = -m (ip(S)gC h), 

hence it must be f^a(hC g) = -fia(C h- g), and the same for B. We can thus 
use Lemma IB. 21 applied to, for each I e A, a scalar product ((■,■)) ^ defined 
in ( 1B.2P using the operators Ck t i,Bi. By Lemma IB. II such a scalar product 
is equivalent to the one 

{hg)j#>\ ■= {h,g) a + V a)l {h ) g). 

a,l — — 

Then Lemma [B.2I states, that 

(h,L {l} h))a,i < -to- 2 (V &l (h,h)+V & i(C o h,C h)) < -ra 2 Va,i(h,h). (6.3) 
where is the generator of the dynamics of the isolated atom I, i.e. 

Lq = £/eA-k{i}- 

The last piece of information we need is given by the following Poincare 
inequality. 

Lemma 6.2. There exist a constant Kp > such that for each a, if f e M^, 
such that Ha(f) = 0, then 

WfWUKp^-DaAfJ)- 

teA 

Proof. By the change of variable introduced in section [0 one can use the 
Poincare inequality for the sphere, for one particle. After that 

l/lli<ll/-^(/)iii+ii^(/)iii 

<K P D &l (f,f) + \\fx ai (f)- f i aha .(f)\\l + \\ (/)! 
iterating the argument yields the result. □ 
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It follows from ( 16. 3ft and Lemma 16.21 that 



i,L h)) = YjiKhi^h))^ < E^' L wM)a,i " E« C o I i'^,Co,i'M>a ) i 



r2 



(6.4) 



rcr II?. II 2 

Accordingly, for each a and h € J^ 1 , such that fi a (h) = 0, we have 

j t {(P*h,Pth}}a = 2{(pth,L pth)) !k 

" l + irp|A|" 
K^n^iP^Plh))^ 

where, in the last line, we have used first f ]6. 3f) and then Lemma [B. II This 
means that 

(P^h,P^h)),<{(h,h)),e-^ 2t 
and, by the equivalence of the norms, there exists C, r' > such that 

llPo^ll^^Ce-^H/ill^i. (6.5) 

□ 



7. Global regularity 

Next, we need regularity of the solution of ( I5.6P also in directions not 
tangent to the energy surfaces, provided g is smooth. Again we work first 
with only one particle. 

7.1. The transversal direction. A natural direction would be given by 
the normal vector to the energy surfaces 



p-d p + VU{q) -d q 



(7.1) 



M 

but, due to the anharmonicity of the potential U, it turns out to be more 
useful to work with the vector field 

Y := ; 1 [ l -p-d P + ~^-q ■ d\ . (7.2) 

Note that, Y is a smooth vector field away from zero. Observe that in the 
harmonic case (U(r) = |r), Y is parallel to the normal vector field Z 3 . 
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The reason why we consider Y is the following. A direct computation 
shows that Y£° = ^yj£°(q,p), hence the vector field is transversal to the 
energy surface. In addition, 



[Y,C ] = 0, 

[Y,L]£ = - 1 -L^/£%q~p)=0. 



(7.3) 



That is, the commutators are vector fields tangent to the constant energy 
surface. Note that since U"(0) > 0, we can write ^^mx = \ + k (\q\ 2 )i f° r a 
smooth function k such that k(0) = 0. An explicit computation yields 

[L, Y] = -A3 7=== Q ■ d P + 7== P ■ g + k1 ]d 9 (7.4) 

Remark 7.1. iVoie i/iai i/ie vector field (£°(q,p))~^ [L, Y] is smooth on all 
M. 2 " , in particular even at zero. This fact will play a crucial role in the 
following. 

7.2. Transversal regularity. The basic idea to prove regularity is to no- 
tice that if Lu = g, then one expects that LYu = [Y, L]u+Yg. Unfortunately, 
we have only L 2 bounds for the right hand side of the above equations, in 
particular we do not know if it belongs to Jif 1 . Hence, a priori, we do not 
know if such an equation has a solution in Jt? 1 . To overcome such a dif- 
ficulty several direct strategies are possible. For example one could try to 
prove a spectral gap in Sobolev spaces of higher regularity or to prove that 
the semigroup maps L 2 functions in Jif 1 functions. Unfortunately, such 
results (even if probably true) are not so easy to prove, in particular the 
related algebra becomes quickly very messy. Due to this state of affair we 
take a bit more indirect route that, without proving explicit bounds, suffices 
to prove the smoothness. To this end it is convenient to work in coordinates 
in which all the energy surfaces can be naturally identified. 

It is then natural to transform equation (14.61) in the following coordinates. 
Let S 3 := {x e M 4 : ||x|| = 1} and M = RxS 3 c RxR 2 xl 2 , M + = (0, oo) xS 3 c 
M. Let : M 4 I A I s {0} -> M[ A] c R5|A| be defined by 




ft/fa), , *^ M ^) 

(7.5) 

The needed properties of this change of variables are detailed in section O 
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The first key observation is that the problem is now regularized at zero 
energy. Indeed, 

ji,k (t, 6, v) = kk ° fc i, v) 

= -j=VV (r,0(rf - r k 9{rie k )i k ) • (r iVi + r kVk ). 

Remark 7.2. Note that j itk extends naturally to a smooth function on M A . 
Indeed if r% < 0, then we can set ji, k (r,£,,r]) = ji, k (r' ,n') where rj = r'-, 
£j = r/j = rj'j for all j + i and ri = -r[, ^ = -^[, r\i - -!][, and the same for 
k. 

Note that for each function / € C°°(Ml A l,R) holds C (f o = (C f) o ^. 
It follows that the equation L$f = g on (R 4 \ {0})l A l is transformed in the 
equation L Q f = g on ^((R 4 n {0})I a I) c M where L = E^C^ + S, and 

f=fo^f-\g = go^-\ 

It is then natural to study the equation on M 

L Ui, k = ji, k . (7.6) 
By the previous discussion the solution of f!4.6p in (R 4 \ {0}) A is given by 

Ui,k = Ui, k \ M A o ^. 

The problem of the transversal smoothness is then reduced to studying 
the smoothness of Ui }k in r (see Lemma [C.ip . 

Lemma 7.3. For each i, k, the functions u i)k e C°°(M A ,R). 

Proof. We can consider Ci , B\ as vector fields on S 3 ^ H Accordingly, we 
can define for each r e R A 

^o(i) = y{° 2 cI + + r(rf4 2 )^]} . 

Then, setting ji tktL (C,,i]) := ji, k (r,£,,r]) and considering the equation on S 3 ^ 

Lo(L)Ui,k,r = ji,k,r (7.7) 

it follows u ik (r,^,r]) = u ik ^,r]). By the previous section we know that, 
for each r, Ui tk . L e C°°(S' 3 I A I,R), thus the differentiability boils down to show 
that the solution of (17. 7p are differentiable with respect to the parameter r. 
Let us fix r and consider the equation^! 

Lo(r)vi, k ,i,r = 9i,k,i, L (7.8) 
2 See ((C21) for the definition of &?. 

3 This is nothing else than the formal derivative of (|7.7p with respect to r; . 
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where gi, k ,r(^v) = ([Y h L ]u i>k ) (r,£,rf) + (Yij i>k )(r,^rf)- Clearly g iA i :L e 
C°°(5 3 I A I,M) for each choice of i,k,l,r, thus u ijJM)E e C°°(S 3 I A I,M). We claim 
that d ri Ui, k = Vi,k,i,r, let us prove it. 
For each, small, h e M A we can write 



L (r + h) 



&i,k,r+h &i,k,r_ ^i,k,l,rh[ 



! 



~ji,k,r+h ji,k,r 9i,k,l,rhl 

- [L (r + h) - L (r)]u i>k , 



An explicit computation shows that the Jif 1 norm of the right hand side 
is bounded by a constant (depending in an unknown manner form r) times 
o(||/i||). The by Proposition 16.11 the differentiability of Ui^, L follows. 

Note that in the above argument the only relevant property of ji : k,r is the 
smoothness on S 3 ^ of itself and of its derivatives with respect to Y[. Since 
a direct computation shows that the same properties are enjoyed by (j^fc^rB 
the Lemma follows by iterating the above argument. □ 



Lemma 7.4. It holds true L^Ui^ € L^ oc (lR 8 , mo). 

Proof. Recall the definition of given by (14.81) . Then all we need to prove 
gisinLfJ 

a* 



is that is in L^ oc (IR 8 , m ). It follows from straightforward calculations 



that the singularities at of are mo integrable. □ 

8. Structure and regularity of 7 2 and a 

Let u = ui t 2 the solution of the Poisson equation (14.61) . Then, by (14. 191) 
and (1D.6I) . we write 

7 2 (ai,a 2 ) = fi ai ,a 2 (ji,2«i,2) 

16ujaia2 f ... v . (8.1) 



/ da (Ci,Vi) do- U 2 ,V2W\/a,C,v)- 

J S 3 xS 3 



Z(a 1 )Z(a 2 ) 

where a is the uniform probability measure on the sphere S 3 and 

(ii,2Mi, 2 )(^ 1 (n,6^i),^" 1 (^,6,?72)) 



Q(r,£,,r)) :-- 



(ii,2Mi,2)(ri,^i,77i,r 2 ,6^2) 
U\p{rl£))u\p{rl$)) 



4 The smoothness on 5 3 ' A ' follows from the smoothness of w^fc, the differentiability 
with respect to Y\ follows by the smoothness of ji,k,r and the fact that Yi'[Yi, Lo]ui^ = 
\Yi>, [Yi, Lo]]ui^ + [Yi, Lo\vi,k,i' ,r, since it is easy to check that, given any vector Z tangent 
to M |A| , [Y V ,Z] is still tangent to M. 
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We have already seen that u satisfies (17. 6j) and is a smooth function on 
M 2 . Hence Q e C°°(M 2 ,R). By Remark El and the structure of L (see 
Lemma IC.lj) , we have 

Gx,2Ui, 2 ) {-n , -6 ,-r/i,r 2 ,C2,^2) = 0i,2«i,2) , 6 , m , r 2 , £ 2 , 7/2) 

and the same for the second coordinate. By the symmetry of the measure 
a it follows then that the integral on the right hand side of (18. ip is an even 
smooth function of y/a[, y/aa, hence a smooth function of c&i, a?,. This shows 
that 7 2 eC°°([0,oo) 2 ,M). 

We are now in the position to prove the relation between a and 7: 

Lemma 8.1. For any nearest neighbor couple {i,k}: 

e u ^(d ai - d ak ) (e-wGV^a*)) = a(a h a k ). (8.2) 

with U (a) = - £7 \ogZ(aj ) . 

Proof. By Lemma ID. II follows 

^W) = e- w a a {eV(^/)} 
provided g does not depend on p. Thus, since m >k is locally bounded, 
a(ai,a k ) = ;u a (VV(ft - q k ){d Pi u i)k - d Pk u i>k )) 
= e~ u d ai {e u Ha(yV(qi - q k )piU itk )} 
- e~ u d ak {e u /j 1 a(VV(q i - q k )p k u i>k )} . 
To continue, notice that for each smooth ip, 

J e -Eio, Y[Z(ai)ip(a)fia(PiVV(qi - q k )u ijk ) 

= E mo ((p(S°)piVV(qi - q k )u ijk ) 

= -E mo (ip(£ )ui tk LnV(qi - q k )) + E mo (ip(£°)p k VV(qi - q k )u i)k ) 
= -E mo ((p(8 °)V(qi - q k )L Ui, k ) + E mo (<p(£°)p k VV '(qi - q k )u iik ) 
= -E mo ((p(£°)V(qi - q k )ji, k ) + E mo ((p(£°)p k VV(qi - q k )u itk ) 
g-Ei^i Y\Z(a l )Lp(a)fXa(p k VV(q i - q k )u ijk ) , 



where we used the antisymmetry of ji jk with respect to p. 
Hence 

Va(PiW(qi - q k )u ijk ) = Ha(p k W(qi - q k )u ijk ) (8.3) 
m -alpmost surely. Accordingly, 

a(a h a k ) = \^e u (d ai - d a J {e u /ia(W (qi - q k ){pi + p k )u iyk )} 



e u (d ai - d ak ) {e u fia (j ik u^ k )} . 
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The result follows since 

Ha{ji,k U i,k) = Ha(Ui,kL U itk ) = a 2 Ha{u i)k Su iik ) 

= ~<y 2 X! ^a((XjU itk ) 2 ) = -^ 2 {a u a k ). 

j=i,k 



□ 



We can rewrite (18. 2p as 

a(a h a k ) = {d ai - d ak )~f 2 (a t ,a k ) + ( -Lr-y - ^7^7 1 7 2 (q», «fc)- (8.4) 

\ Z[ai) Z(a k ) J 

Since ~ a -1 , the regularity of a follows from the one of 7 2 , if we can 
prove that 7 2 ~ 

Lemma 8.2. There exists G e C°°([0, oo) 2 ,IR) such that, for each 01,02 > 



7 2 (ai,a 2 ) = a 1 a 2 G(ai,a 2 ) > 0, 
and furthermore a(0, a) > /or all a>0. 

Proof. Observe that fixing the energy of the first particle £i(qi,Pi) = 0, it 
implies that qi = pi = 0. So defining u(q,p) = 1*1,2(0, 0, q,p), it solves on M. 2 
the equation Lu = ], where L is the generator of the dynamics of a single 
isolated atom and j = -\p ■ VV(q) = L*V . 
By the smoothness of 7 2 it follows that 

7 2 (0, a) = -n a Qu) = {i a (L*V ■ u) = fi a (VLu) = fi a (Vj) = 

due to the symmetry of /i a with respect to the transformation p -»■ -p. The 
structure of 7 2 follows then by the symmetry and smoothness of 7 2 . The 
positivity follows by 

7 2 = -H ai ,a 2 (L U-u) = (T 2 /i aiia2 ((XiM) 2 + (X 2 u) 2 ). 

By dH3J and flTXo]) 

a(a,0) =2aG(a,0) > 0. 

□ 

Appendix A. Commutators 

This appendix collects some formulae concerning commutators for the 
unperturbed system (e = 0), we use the notation of Proposition 16.11 In 
section [5] we have already computed: 

[C , B] = (Jp-d q + JVU(q)-d p ) =: C x . (A.l) 
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(A.2) 



[d, B] = AU (q)Jq ■ d q - A{U (q)\q\ 2 + U}Jp-d p 
=:2C 2 + C(g)Co=:2C 2 + J R 2 , 
An explicit computation shows that 

C2 __PjJvu B + r^u Ci _mL Co . (A . 3) 

'jp 'P *p 

The above formula shows that C 2 , for small p and large q, is not bounded 
by B,Co,Ci and, as we will see in the following, this forces us to compute 

R 3 = [C 2 , B] = -{JVU(q) ■ D 2 U(q)d p + JD 2 U(q)p ■ d q ) 

= -{4U"{p, q)Jq + 2U'jp}d q - 4(U') 2 Jqd p . 

By using the orthonormal base Z ,Zi,Z 2 , defined in (|5.5p . we have that 

2 
3=0 

= Af- 2 (R 3 , B)B+ AA- 2 (i? 3 , C x )Ci + AT 2 (i? 3 , (C 2 - C ))(C 2 - C ) 
= 4Af- 2 U"(q,p)(Jp, q)B-M- 2 \4u" (q, p) 2 + 2U'p 2 + 8(17') 3 g 2 ] C x 

-8Af- 2 u"u'(q,p)q\C 2 -C ), 

and an explicit calculation shows that 

||J?s|| 2 = AT 2 (\(R 3 , B)\ 2 + \(R 3 , C 1 )\ 2 + \(R 3 , (C 2 - C ))\ 2 ) < K (A.4) 

To conclude the first order analysis we need to compute some more com- 
mutators 

([C h C 1 ]) = -B 

([C ,C 2 ]) = (A.5) 
([Ci,C 2 ])=: pC 2 -(5B, 

with p and /3 also bounded. 

We also need some second order commutators: 

[B, Cq] = CqBCq - CqB - C\Cq = -2C\Cq + B 

[Ci, Cq] = C B + BC = 2C B - C\. 
Appendix B. Hypocoercivity estimates 



(A.6) 



This appendix contains the core of the hypocoercivity argument. For our 
purposes it turns out to be more convenient to set it in an abstract setting. 

Let VP be an Hilbert space and Co,B be closed operators satisfying the 
relations (1A.1I) - (1A.6I) . Assume that B,CkCo, BCq have all a common 
core D c . In addition, assume that for each h,g e D c and Z e {Cq,B}, 

(h,Zg) = -(Zh,g), (B.l) 
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where (•,•) is the scalar product of "H . Remark that the various constant 
that will appear in the results (t,K etc.) do not depend on the scalar 
product (h,g). 

We are interested in obtaining estimates in terms of the following gener- 
alized Sobolev norm: 

2 

\\h\\%i ■■= m 2 +T,t c i> h t 2+ t Bh t 2 - 

k=0 

To this end it turns out to be convenient to define the bilinear form: 

2 

((h,g)) -={h,g) + Y a k( c kh,C k g) + a 3 {Bh,Bg) 

k=o (B.2J 

- b (C h, C x g) - boiCth, C g) - ^{dh, C 2 g) - h{C 2 h, C x g) 
where a k > 0,b k > will be chosen shortly (see flB.6j) and (IB.7|) ) . 
Lemma B.l. If for 5 e (0, 1) 

a>obj + a 2 b\ < a aia 2 (l - 5) 2 , (B.3) 
then the quadratic form is positive definite and 

«/ i ,/ i ))>(||/ i || 2 + 5^a fe ||C fe / i || 2 + a3|| J B/ i || 2 ) = K||/ i ||^ 1 . 

fc=0 

with k = min{5ao, Sax, 5a 2 , 0,3, 1}. 

Proof. We have for any cto, ot-x > 

2 1 

((h, h)) - a 3 \\Bh\\ 2 - \\h\ 2 > Y, a k \\C k h\\ 2 - Y (b k a k \\C k h\\ 2 + ha^WC^xhf) 

fc=0 k=0 

= (a - b ao) \\C h\\ 2 + (ax - fro^o 1 - ^i^i 1 ) |Ci^| 2 + ( a 2 - bx&x) \\C 2 h\\ 2 
2 

=5 Y a k \\C k h\\ 2 + (a (l - 5) - b a ) \\C h\\ 2 + (a 2 (l - 5) - b x a x ) \\C 2 h\\ 2 

fc=0 

+ (ai(l - 5) - feotto 1 ~ ^i"! 1 ) |Ci^| 2 
Then choosing 

a (l-5) a 2 (l-8) 
a = , ax = 1 

b bx 

the Lemma follows immediately by condition flB.3j) . □ 

On the other hand Schwartz inequality implies that there exists K > 
such that 

{{h, h)) < K(\\h\\ 2 + Y \\C k h\\ 2 + \\Bh\\ 2 ). (B.4) 

fc=0 

Let % l ■= {h e 7-L° : h} < 00}, clearly it is an Hilbert space, equivalent 
to Jf? 1 , with scalar product ((•,•)). 
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Lemma B.2. Given Cq,B as described at the beginning of the section, there 
exists t > such that, for each a e (0, 1), 

(h, Lh) < -to- 2 [ £ \\C k hf + £ ||C fe C /*| 2 + imi 2 + liWl 2 } . (B.5) 

U=0 fc=0 J 

Proof. This is a proof by boring computations. Let us start: 
{{h, Lh} =a 2 (h, C 2 h) + { £ a fc [a 2 <C fe C 2 /i, C fe /i) + (C k Bh, C k h)] 

lA.k lB,k 

+ a 2 a 3 (BC 2 h, Bh) - £ a 2 6 fc ((C fc /i, C fe+1 C 2 /i> + <C fe C 2 /i, <7 fc+1 fc)) 



fc=0 



Ib.b II A) 

X 6, ({C k h, C k+l Bh) + (C k Bh, C k+1 h}) }. 



k=0 



Hb,] 



Now we must look at all the terms one by one, we will use systematically 

(h,C 2 h) = -(C h,C h). 

(C*h,C h) = -(C 2 h 1 C 2 h). (I Ai0 ) 

(CiCg/i, C^/i) = (QACo/i, Cift) + (SC /i, C x h) 

= -(dCoh, Codh) + (C Bh, dh) - {dh, dh) (Ia,i) 
= -{dC h, dC Q h) + (CiC /i, - (S/i, C Ci/i) 

~{dh,dh) 
= -(dC h, ddh) + (Bh, Bh) - (dh, dh). 

(dC 2 h,dh) = -(ddKddh). (I A)2 ) 

(C Q Bh,C h) = (dh,C h). (I B ,o) 

(dBh,dh) = 2(dh,dh) + (R 2 h,dh). (I B ,i) 



(dBh,dh) = (R 3 h,dh). 



2(» 
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(BC 2 h, Bh) = -(BC h, C Bh) - (CxCoh, Bh) 

= -{BC h, BC h) - {BC h, C x h) - (ddK Bh) (I BjB ) 
= -(BC h, BC h) + (Bh, Codh) + (dh, dh) - (ddh, Bh) 
= -{BC h, BC h) - {Bh, Bh) + (dh, dh). 



(dC 2 h, dh) + (C h, CtClh) = -2(C 2 h, dC h) - (C h, dh). (II A , ) 



(dC 2 h, dh) + {dh, dC 2 h) = -2(dC Q h, ddh) - {dh, dh). (II A ,i) 



{dBh,dh) + {dh,dBh) = {dh,dh)+2{dh,dh) + {dh,R 2 h). (II B) o) 



{dBh, dh) + {dh, dBh) = 2{dh, dh) + (R 2 h, dh) + {dh, R 3 h) . (II B ,i) 
Finally, we can put all the terms together, obtaining 

(h, Lh)) = - a 2 \\dh\\ 2 - [b + a 2 ( ai - a 3 )] \\dhf - 2b, \\dh\\ 2 

- a 2 (a 3 - a^WBhf + (a + a 2 b ) (dh, dh) + a x (R 2 h, dh) 

- b (R 2 h, dh) - h(R 2 h, dh) + (2ai + a 2 h)(dh, dh) 

- 2b (dh, dh) + a 2 (R 3 h, dh) - h(R 3 h, dh) 

- a 2 ( £ a k \\ddh\\ 2 + a 3 \\Bdh\\ 2 - 2 £ b k (ddK C k+1 C h)\ 

\k=0 k=0 ) 

v ' 

III 

By the same argument used in the proof of Lemma IB.lt if for < 5 < 1 we 
have 

a bl + a 2 bl < a ctia 2 (l - 5) 2 , (B.6) 

then III is bounded by -a 2 S (£* a k \\C k C Q h\\ 2 + a 3 \\Bdh\\ 2 ). 

Recalling flAJ2]) and (EOj) . we have \\R 2 h\\ 2 < K 2 \\dh\\ 2 , \\R 3 h\\ 2 < K 3 {\\Bh\\ 2 + 
+ \\Cih\\ 2 + \\dh\\ 2 }. Applying Schwarz, for each a%, a 2 , a 3 , > 0, 
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the first three lines of the above equation are bounded by 

- jcr 2 - ^«i[a + a 2 b + K 2 a x ] - K 2 (b + b x ) - a 3 b - ^K 3 a 2 - -hK 3 a^ \\C h 

f 2 a + a 2 b + a 1 2a x + a 2 h K 3 a 2 {K 3 a\ + ) .. 
-r°- a(fl3 - ai) 2^ 2^ 2 2^ 

- |26i - -6i - ia 2 (2a! + a 2 6i) - a 3 \ - ^(K 3 + l)a 2 - i&^c^J |C 2 ^| 2 

- ja 2 (a 3 - ai) - ^X 3 a 2 - ^K 3 b x a^ \\Bh\\ 2 

An inspection of the above expression shows that with the choices 

a = a 2 v 7 ; a x = a 2 v 16 ; a 2 = a 2 v 18 ; a 3 = v 13 ; b = a 2 v 12 ; bi = a 2 v 17 (B.7) 
and 

«i = v' 6 ; a 2 = v 2 ; a 3 = v~ 5 ; a 4 = v 

implies that, by choosing v small enough, there exist r > such that, for 
each a e (0, 1), 

(h, Lh) < -ra 2 I £ \\C k h\\ 2 + \\Bh\\ 2 + £ \\C k C h\\ 2 + \\BC h\\ 2 \ . (B.8) 

U=0 k=0 J 

Observe that (IB. 61) is satisfied by this choice. □ 

Appendix C. A coordinate change 

We study a change of coordinates in the case on one particle, for many 
particles one simply considers the product. 

Let S 3 := {x € R 4 : ||a;|| = 1} and M = R x S 3 c R x R 2 x R 2 , M + = 
(0, oo ) x S 3 c M. Let \I> : R 4 \ {0} -» M + c R5 be defined by 

/ 



Remember that, by hypotheses, = U(q 2 ) and {7 is a strictly increasing 
function. We can then extend U to a smooth increasing function on R such 

that 77(2) > if z > OH It follows that p(z) := U~ (z) is a well defined 
smooth function on R such that p(0) = 0. 

One can readily check that the inverse ^J^ 1 ■ M + -*■ R 4 is given by 



(C.l) 



3 Clearly the extension is arbitrary, but this is irrelevant in the following. 
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For the following it is convenient to introduce the function 9(z) = \J p(z)/z, 

2^0. Notice that 9 is smooth on K provided we set 9(0) = \Jl/U (0). 

Next we transport the vector fields on M by the usual formula ^>*Z = 
(D^/Z) o^ _1 . The following lemma follows by the computation of that 
can be found in appendix ID1 

Lemma C.l. With the above notations we have 
C := tf,C = Jr]-d v 

-V2f7'(p(r 2 e 2 ))^r 2 e 2 )^ 

Y := V*Y = d r . 
Define the function 

T(z) :=1-Tj'(p(z))9(zf 

and notice that it is continuous in and T(0) = 0. It is then natural to 
write 

U[r * > (C.2) 

Notice that 9(r 2 ^ 2 ) and r(r 2 £ 2 ) are smooth function of £ and r, conse- 
quently B is a smooth vector field. In addition, both B° and B 1 are tangent 
to the surfaces {r} x 5 3 c M. 

Remark C.2. The vector fields Cq,B,Y are defined only on M + but admit 
a smooth canonical extension to all M. We will use the same notation to 
designate such an extension. 

It follows that there exists smooth vector fields C\,Ci on M such that 
Q = #*C* on M + . 



Appendix D. Microcanonical measure 

We collect here some properties about the microcanonical measure p, a 
(and its product version Ha). 

Let us first recall the definition and some formulas. Microcanonical mea- 
sure is defined as the conditional measure on the energy shell S a = {(q,p) £ 
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R 4 : |p| 2 /2 + U(q) = a}. This means that for every continuous functions 
<p : R + -+ R and / on R 4 we have 

f 4>(£o(p,q))f(p,q)dm = f 4>(a)^ a (f)e- a Z(a)da (D.l) 

J JR+ 

Standard formulas give that 

2(a) = f ( P 2 + U\ q Yy ll2 da^(P,q) (D.2) 

and 

Hatf) = Z{a)- 1 f f{p,q){p 2 + U>{ q yy ll2 d^ a {P,q) (D.3) 

where a^ a is the Lebesgue measure on the energy shell. 

We need an integration by part formula for the microcanonical measure 
/i a . This is provided by the following proposition. 

Proposition D.l. For all f,g continuous locally bounded functions on R 4 ; 
such that g does not depend on p k (d Pk g = 0), and f is differentiate in pk, 
then 

Z' ' ( 1 

Va(gdpJ) = (d ah Ha)(gp k f) + — k , Ha(gPkf), m 0>r a.s., a k > 0. (D.4) 

Z{a k ) 

Proof. Because of the product structure of we can just consider only one 
site and we drop the index A; of it. 

Let us recall that from the definition of fJL a follows, for each smooth 
bounded functions /, g and <p : R+ -*■ R supported away from zero, 

f dpdqe- H ^<p(H )f(p,q) = f da e' a ip{a)Z{a)^i a {f) 

JR 4 JR+ 

Thus 

/ da e- a (f(a)Z(a)^ a (gd p f) = dp dq e- H °ip(H )g(q)d p f(p, q) 

JR + JR 4 

= f dpdqe- H °{-cp>(H ) + tp(H )}g(q)pf(p iq ) 

JR 4 

= / dae- a {-v'(a) + ip(a)}Z(a)ii a (gpf) 

JR+ 

= / dae- a (p(a)d a {Z{a)n a {gpf)) 

JR+ 

= / dae- a ip(a)Z(a)d a ({i a (gpf)) + / da e~ a ip(a) Z(a)fi a (gpf). 
Jr+ Jr+ Z\a) 

It follows the relation (ID M for any bounded smooth /, g. The result follows 
by approximations. □ 
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Formula ( ID. 3}) is difficult to be used directly, but exploiting the symmetry 
and the convexity of the potential, it is possible to write this microcanonical 
expectation as as integral on the 3 dimensional sphere of radius 1 with 
respect the corresponding uniform measure. In fact the strict convexity of 
U makes the energy shell very close to a sphere for small energy a, and p a 
close to the uniform measure on this sphere. We want to study this more 
precisely. 

Recall the change of coordinates ^ introduced in Section by ( I7.5P 
and its inverse (jC.ip . Recall also the notation 9(z) = \J p(z)/z, and that 

6(0) - VW(0)- 
Lemma D.2. 

Z(a) = 4u 4 a f^Uipiae))}' 1 daU, V ), (D.5) 

where we have used polar coordinates in four dimensions, a is the uniform 
probability measure on S 3 , the unit four dimensional ball, and u; 4 is its 
volume. Furthermore 

Pa(f) = ^§)J S J° p'ipiV^e))] 1 da&ri). (D.6) 

Proof. Instead of computing with differential forms it turns out to be more 
efficient to use the following trick. 



Consider the change of variables \I/ : (R 4 \ {0}) 



I 5 defined by 



s q 



sp 



\q\y/^ + U(q) \/p 2 + 2[/(<?) 
Note that, ^ is invertible and, setting g(£,fj) = \/f 2 + rj 2 , 



(D.7) 



^\Jp(ag 2 £ 2 )-^,£ 1 V2afj,Q^ 

Then, given any two compact support functions g £ I 
we can write 



/€C°(I 



g(s)f(p,q)dqdpds= I \gf \det(D^/)\ oty 1 (a,£,fj)dad£dr] 

$5 Jro l 1 J 

= 4w 4 J ^ s g(s) [^/|det(D§)p 1 o^- 1 (a,s^,sr])s 3 dadsda(^,r]) 



In Lemma ID. 3 1 we compute the determinant of 



{ da da da \ 

dq dp ds 

d£ d£ d£ 

dq dp dp 

dr\ dr\ dij , 

\ dq dp dp I 



(D.8) 
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for £ 2 + rf = 1, obtaining 



det(Mo*-M=^ 2 » 

V ' rt 



(D.9) 



Thus, if we take a sequence of g n that converges weakly to the delta 
function on one, we have the formula 

J(p,q)dqdp = 4u 4 / f o^-^a,^!]) - dada(£,r)). 

Accordingly, for each g e L 1 (R + ) we can write 

f(p,q)g(£°)dqdp = 4u 4 / g(a)f o qHfo^) ° dadcr(£,r]) 

On the other hand by the definition of the microcanonical measure /i a : 
/ f(p,q)g(£°)dqdp= / g(a)Z(a)fj, a (f)da, 

JR 4 JM.+ 

The above, by the arbitrariness of g, yields the following formula for the 
microcanonical expectation: 

fjt a (f) = ^aZ(a)- 1 f f o ^-\a,C,v) [^'(pK 2 ))]^ M^V)- P-IO) 

Putting / = 1 in fID.lOp implies the formula for Z(a). □ 
Lemma D.3. Proof of equation (1D.9|) . 
Proof. An explicit computation of the derivative yields 

We want to compute the determinant for (£,,rj) e S 3 , i.e. g = 1. If we 
multiply the first row by (2a) _1 £i and we sum it to the second, then by 
(2a) _1 ^ 2 and sum it to the third, by (2a) -1 ?7i and sum to the fourth and 
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finally by (2a) 1 rj2 and sum it to the last row, we have 

/ 2£/'0« 2 ))#« 2 )v/< V2a~r] 0\ 



det o |r-i) = det 







o e 



l(2a)~2 



/ 2U'(p(ae))e(ae)^ 



det 





l(2a)~ 



-2ary 2 \ 

e 



= — det 
2a 



' 2£/'(p(a£ 2 ))#« 2 )v / a£ -2ar? 2> 



2a 



det 



-2a?7 2 \ 



Prom the above the Lemma follows]^ □ 

Corollary D.4. Let f be a continuous function of(q,p), then the following 
formula holds 

Mq,p))= ffo^a^da^ + O^a) (D.ll) 
j s 3 

where a is the uniform probability measure on S 3 , and 0(a) is a smooth 
function of order a as a -> 0. 



Notice that because of ( ID. 101) . since < c < U' < cr 1 < +oo, microcanonical 
measure is uniformly equivalent to the uniform measure on the unit sphere, 
and, for any a > 0, and any positive /: 



,-2 



f f o (a, e, 7/) V)<^a(f)< c 2 f f o (a, e, /?) d<x(£, 77) 

(D.12) 



6 Just remember that det(l -bv ®v) = 1- bv 2 since i> and any vector perpendicular to 
v are eigenvectors of eigenvalue 1 - bv 2 and 1 respectively. 
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Appendix E. The Gaussian case 

Just to give a concrete idea of what we are doing and to provide some 
concrete intuition, here we discuss the case in which we have just two parti- 
cles and both U and V are quadratic. The point being that such a case can 
be solved explicitly and hence provides a guidance for what can be expected 
in the general case. 

Let us consider and Hamiltonian system with four degree of freedom 
(?i, Pi, 92,^2) e given by 

1 1 
H £ ■= ^{\\Pi\\ 2 + \\P2W 2 } + ^{\\qi\ 2 + IMI 2 } + e\ qi -q 2 \ 2 

plus random forces that conserves the kinetic energy (that is independent 
diffusions on the circles ||pi|| 2 = cost). To this end consider the vector fields 

Xi -P^idp^-pi^dp^. 

The generator is thus given by 

L £ :={H £ ,-} + a 2 Y j Xf = A e + a 2 S 

i=i 

The energies of the two particles are £i = |||pi|| 2 + \\qi\ 2 + §||<?i - (fell 2 and 

d t £i = e(q 2 -q 1 ,p 1 + P2)- ej 
d t £2 = ~ej 

gives the current. A direct computation shows that, setting 
u '= \{ II ?2 1| 2 - || 9i || 2 } - o-- 2 q 2 ■ pi + (J- 2 q x ■ P2 

holds true 

L £ u = j + Ea- 2 {\\ qi \\ 2 -\\q 2 \\ 2 }=:j. (E.l) 

Accordingly, if we rescale the time by e~ 2 and we look at the random vari- 
ables £i, £ (t) '■= £i(s~ 2 t). 

M*)-M0) = * f 2 \Leu){s)ds-e 2 a- 2 H ^ {\\ qi {s)\ 2 - \\q 2 {s)f)ds 
Jo Jo 

= e (u(e- 2 t) - u(0)) + eMl H - a~ 2 f\^Ar) ~ ^, e (r))dr 

-^ 2 £ 2< [(ll9i(^)ll 2 -bi(^)ll 2 )-(lk2( s )ll 2 -b 2 (^)ll 2 )]^^(^ 3 ) 

(E.2) 

where the martingale M" has quadratic variation given by 
< M u > t = a 2 £ [(X lU (s)) 2 + (X2u(s)) 2 ] ds 
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We first show that the average of the last term on the RHS of ( IE. 21) tends 
to 0. Observe that 

(Ikill 2 -bill 2 )-(lk2|| 2 -b2|| 2 ) 

= -L B (pi -qi -p 2 -Q2 + ^(1 gill 2 - Ml 2 )) - e(M| 2 - Ml 2 ) 

Calling v = pi ■ qi - P2 ■ q2 + ^(\\qi\\ 2 ~ \\q2W 2 ) , the last term on the RHS of ( ]E.2j) 
can be rewritten as 

e 2 (v(e-H) - v(0)) + e 2 Ml_ H + e 3 £ ^ (\\ qi (s)\\ 2 - \\q 2 (s) \\ 2 )ds 

It is easy to show that the average goes to as e -> 0. 

It remains to compute the limit of the martingale eM™ 2( . To this purpose 
one has to compute the limit of its quadratic variation 

[(<?2,lPl,2 - ?2,2Pl,l) 2 + (<?1,2P2,1 - qi,iq2,2) 2 ] ds 

After some explicit calculation, one can find 4th order polynomials ip,<f),x 
such that 

(?2,lPl,2 - <?2,2Pl,l) 2 + (<?1,2£>2,1 - <?1,1?2,2) 2 = 2£i£ 2 + + S(f) + S\ 

The term L £ ifj + ex does not contribute. For the term S<j), by using Schwartz 
inequality: 

~fo dT I S( l )F r d ^ = fdr^ J (XiMXiFjdfte 

<[£^F T )dr)j 1/2 (£^f {X4) 2 F T d^ E dr\ 



£ 1 /2 / t ~ \ 

<^*X>(F T )dr) f jf £ (^0) 4 rf/i £ ] 2 rfrj < CS/t 

We have obtained that £i >£ ,£2,e converge to the (degenerate) diffusion on 
Rl generated by 

C = 2a- 2 (d £l - d e ,)EMd £l - d £2 ) (E.3) 
corresponding to the stochastic differential equation: 

dSi = a' 1 xllExEi dw t - 2a' 2 (S 1 -£ 2 )dt= -d£ 2 (E.4) 
with wt standard Wiener process. 
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